Crop Phenology Estimation Using a Multitemporal Model and a Kalman Filtering Strategy

In this letter, a new approach for crop phenology estimation with remote sensing is presented. The proposed methodology is aimed to exploit tools from a dynamical system context. From a temporal sequence of images, a geometrical model is derived, which allows us to translate this temporal domain into the estimation problem. The evolution model in state space is obtained through dimensional reduction by a principal component analysis, defining the state variables, of the observations. Then, estimation is achieved by combining the generated model with actual samples in an optimal way using a Kalman filter. As a proof of concept, an example with results obtained with this approach over rice fields by exploiting stacks of TerraSAR-X dual polarization images is shown.This project was supported in part by the Spanish Ministry of Economy and Competitiveness (MINECO) and in part by EU FEDER under Project TEC2011-28201-C02-02

Dynamical Approach for Real-Time Monitoring of Agricultural Crops

In this paper, a novel approach for exploiting multitemporal remote sensing data focused on real-time monitoring of agricultural crops is presented. The methodology is defined in a dynamical system context using state-space techniques, which enables the possibility of merging past temporal information with an update for each new acquisition. The dynamic system context allows us to exploit classical tools in this domain to perform the estimation of relevant variables. A general methodology is proposed, and a particular instance is defined in this study based on polarimetric radar data to track the phenological stages of a set of crops. A model generation from empirical data through principal component analysis is presented, and an extended Kalman filter is adapted to perform phenological stage estimation. Results employing quad-pol Radarsat-2 data over three different cereals are analyzed. The potential of this methodology to retrieve vegetation variables in real time is shown.This work was supported in part by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER under Project TEC2011-28201-C02-02 and in part by the Generalitat Valenciana under Project ACOMP/2014/136

A Simple RVoG Test for PolInSAR Data

In this paper, we present a simple algorithm for assessing the validity of the RVoG model for PolInSAR-based inversion techniques. This approach makes use of two important features characterizing a homogeneous random volume over a ground surface, i.e., the independence on polarization states of wave propagation through the volume and the structure of the polarimetric interferometric coherency matrix. These two features have led to two different methods proposed in the literature for retrieving the topographic phase within natural covers, i.e., the well-known line fitting procedure and the observation of the (1, 2) element of the polarimetric interferometric coherency matrix. We show that differences between outputs from both approaches can be interpreted in terms of the PolInSAR modeling based on the Freeman-Durden concept, and this leads to the definition of a RVoG/non-RVoG test. The algorithm is tested with both indoor and airborne data over agricultural and tropical forest areas.This work was supported by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER under Project TEC2011-28201-C02-02

Influence of Incidence Angle on the Coherent Copolar Polarimetric Response of Rice at X-Band

The coherent nature of the acquisition by TerraSAR-X of both copolar channels (HH and VV) enables the generation of many different polarimetric observables with physical interpretation, as have recently been used for monitoring rice fields. In this letter, the influence of incidence angle upon these polarimetric observables is analyzed by comparing three stacks of images that were acquired simultaneously at different incidence angles (22°, 30°, and 40°) during a whole cultivation campaign. We show that the response of observables related to dominance (entropy, ratios of components) and type of scattering mechanisms (alpha angles) is not greatly influenced by incidence angle at some stages: early and advanced vegetative phases, and maturation. Moreover, the acquisition geometry drives the sensitivity to the presence of the initial stems and tillers, being detected earlier at shallower angles. This analysis is a necessary step before studying potential methodologies for combining different orbits and beams for reducing the time between acquisitions for monitoring purposes.This work was supported by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER under Project TEC2011-28201-C02-02, and by the Generalitat Valenciana under Project ACOMP/2014/136

Particle Filter Approach for Real-Time Estimation of Crop Phenological States Using Time Series of NDVI Images

Knowing the current phenological state of an agricultural crop is a powerful tool for precision farming applications. In the past, it has been estimated with remote sensing data by exploiting time series of Normalised Difference Vegetation Index (NDVI), but always at the end of the campaign and only providing results for some key states. In this work, a new dynamical framework is proposed to provide real-time estimates in a continuous range of states, for which NDVI images are combined with a prediction model in an optimal way using a particle filter. The methodology is tested over a set of 8 to 13 rice parcels during 2008–2013, achieving a high determination factor R2=0.93 ( n=379 ) for the complete phenological range. This method is also used to predict the end of season date, obtaining a high accuracy with an anticipation of around 40–60 days. Among the key advantages of this approach, phenology is estimated each time a new observation is available, hence enabling the potential detection of anomalies in real-time during the cultivation. In addition, the estimation procedure is robust in the case of noisy observations, and it is not limited to a few phenological stages.This work is supported by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER under Projects TEC2011-28201-C02-02 and TIN2014-55413-C2-2-P

Polarimetric Response of Rice Fields at C-Band: Analysis and Phenology Retrieval

A set of ten RADARSAT-2 images acquired in fully polarimetric mode over a test site with rice fields in Seville, Spain, has been analyzed to extract the main features of the C-band radar backscatter as a function of rice phenology. After observing the evolutions versus phenology of different polarimetric observables and explaining their behavior in terms of scattering mechanisms present in the scene, a simple retrieval approach has been proposed. This algorithm is based on three polarimetric observables and provides estimates from a set of four relevant intervals of phenological stages. The validation against ground data, carried out at parcel level for a set of six stands and up to nine dates per stand, provides a 96% rate of coincidence. Moreover, an equivalent compact-pol retrieval algorithm has been also proposed and validated, providing the same performance at parcel level. In all cases, the inversion is carried out by exploiting a single satellite acquisition, without any other auxiliary information.This work was supported by the Spanish Ministry of Economy and Competitiveness (MINECO) and European Union FEDER under Project TEC2011-28201-C02-02

Estimation of Key Dates and Stages in Rice Crops Using Dual-Polarization SAR Time Series and a Particle Filtering Approach

Information of crop phenology is essential for evaluating crop productivity. In a previous work, we determined phenological stages with remote sensing data using a dynamic system framework and an extended Kalman filter (EKF) approach. In this paper, we demonstrate that the particle filter is a more reliable method to infer any phenological stage compared to the EKF. The improvements achieved with this approach are discussed. In addition, this methodology enables the estimation of key cultivation dates, thus providing a practical product for many applications. The dates of some important stages, as the sowing date and the day when the crop reaches the panicle initiation stage, have been chosen to show the potential of this technique.This work was supported in part by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER under Project TEC2011-28201-C02-02, and in part by the Generalitat Valenciana under Project ACOMP/2014/136. All SAR images have been provided by DLR in the framework of projects LAN0021 and LAN0234 of the prelaunch AO of TerraSAR-X

Retrieval of phenological stages of onion fields during the first year of growth by means of C-band polarimetric SAR measurements

The phenological stages of onion fields in the first year of growth are estimated using polarimetric observables and single-polarization intensity channels. Experiments are undertaken on a time series of RADARSAT-2 C-band full-polarimetric synthetic aperture radar (SAR) images collected in 2009 over the Barrax region, Spain, where ground truth information about onion growth stages is provided by the European Space Agency (ESA)-funded agricultural bio/geophysical retrieval from frequent repeat pass SAR and optical imaging (AgriSAR) field campaign conducted in that area. The experimental results demonstrate that polarimetric entropy or copolar coherence when used jointly with the cross-polarized intensity allows unambiguously distinguishing three phenological intervals.This work was partially supported by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER, under Project TEC2011-28201-C02-02

Retrieval of vegetation height in rice fields using polarimetric SAR interferometry with TanDEM-X data

This work presents for the first time a demonstration with satellite data of polarimetric SAR interferometry (PolInSAR) applied to the retrieval of vegetation height in rice fields. Three series of dual-pol interferometric SAR data acquired with large baselines (2–3 km) by the TanDEM-X system during its science phase (April–September 2015) are exploited. A novel inversion algorithm especially suited for rice fields cultivated in flooded soil is proposed and evaluated. The validation is carried out over three test sites located in geographically different areas: Sevilla (SW Spain), Valencia (E Spain), and Ipsala (W Turkey), in which different rice types are present. Results are obtained during the whole growth cycle and demonstrate that PolInSAR is useful to produce accurate height estimates (RMSE 10–20 cm) when plants are tall enough (taller than 25–40 cm), without relying on external reference information.This work has been supported by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER under project TIN2014-55413-C2-2-P. The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7/2007–2013) under grant agreement 606983, and the Land-SAF (the EUMETSAT Network of Satellite Application Facilities) project. The in-situ measurements in the Ipsala site were conducted with the funding of The Scientific and Technological Research Council of Turkey (TUBITAK, Project No.: 113Y446)

Sentinel-1 InSAR coherence for land cover mapping: a comparison of multiple feature-based classifiers

This article investigates and demonstrates the suitability of the Sentinel-1 interferometric coherence for land cover and vegetation mapping. In addition, this study analyzes the performance of this feature along with polarization and intensity products according to different classification strategies and algorithms. Seven different classification workflows were evaluated, covering pixel- and object-based analyses, unsupervised and supervised classification, different machine-learning classifiers, and the various effects of distinct input features in the SAR domain—interferometric coherence, backscattered intensities, and polarization. All classifications followed the Corine land cover nomenclature. Three different study areas in Europe were selected during 2015 and 2016 campaigns to maximize diversity of land cover. Overall accuracies (OA), ranging from 70% to 90%, were achieved depending on the study area and methodology, considering between 9 and 15 classes. The best results were achieved in the rather flat area of Doñana wetlands National Park in Spain (OA 90%), but even the challenging alpine terrain around the city of Merano in northern Italy (OA 77%) obtained promising results. The overall potential of Sentinel-1 interferometric coherence for land cover mapping was evaluated as very good. In all cases, coherence-based results provided higher accuracies than intensity-based strategies, considering 12 days of temporal sampling of the Sentinel-1 A stack. Both coherence and intensity prove to be complementary observables, increasing the overall accuracies in a combined strategy. The accuracy is expected to increase when Sentinel-1 A/B stacks, i.e., six-day sampling, are considered.Peer ReviewedPostprint (published version